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ABSTRACT 

Accurate measurements of the cosmic microwave background (CMB) anisotropics with 
an angular resolution of a few arcminutes can be used to determine fundamental cos- 
mological parameters such as the densities of baryons, cold and hot dark matter, and 
certain combinations of the cosmological constant and the curvature of the Universe 
to percent-level precision. Assuming the true theory is a variant of inflationary cold 
dark matter cosmologies, we calculate the accuracy with which cosmological parame- 
ters can be determined by the next generation of CMB satellites, MAP and Planck. 
We pay special attention to: (a) the accuracy of the computed derivatives of the CMB 
power spectrum C^; (b) the number and choices of parameters; (c) the inclusion of 
prior knowledge of the values of various parameters. 
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1 INTRODUCTION 

The detection of primordial anisotropies in the microwave 
background radiation by the COBE satellite (Smoot et al. 
1992) has had an enormous impact on cosmology (see White, 
Scott and Silk 1994 and Bond 1996 for reviews). How- 
ever, the relatively poor angular resolution of COBE/DMR 
{Ofwhm ~ 7°) limits the amount of information that can 
be extracted from the CMB. From the 4 year COBE maps 
K^ i: (Bennett et al. 1996a), the overall amplitude of the CMB 
» I ■ power spectrum for a given spectral shape has been deter- 
C^ ■ mined to an accuracy of 7% and a power law index char- 
acterizing the shape to ±0.24 (Bond 1996). Constraints on 
other parameters such as the spatial curvature and the cos- 
mological constant A are weak. 

It has long been known (Sunyaev and Zeldovich 1970) 

that at angular resolutions smaller than ~ l°fio i^^^ ^^~ 
gle subtended by the Hubble radius at the time of recom- 
bination) the CMB power spectrum will depend on e.g., 
the sound speed of the baryon-photon fluid, and hence on 
a number of fundamental cosmological parameters, such as 
the densities of baryons, cold and hot dark matter, and the 
spatial curvature of the Universe. In adiabatic models, the 
acoustic motions of the matter radiation fluid lead to a char- 
acteristic series of 'Doppler peaks' in the CMB power spec- 
trum which have been investigated in considerable detail 
numerically and semi-analytically (e.g.. Bond 1996, Hu et 
al. 1997). Similar behaviour is expected qualitatively in de- 
fect (isocurvature) theories, though the pattern of Doppler 



peaks is expected to be less distinct and has not yet been 
calculated to high precision (e.g., Turok 1996). We therefore 
restrict the discussion in this paper to purely adiabatic per- 
turbations obeying Gaussian statistics, as expected in most 
inflationary models of the early universe (e.g., Linde 1990). 
The anisotropies in such models can be computed to high 
accuracy which, as we will show in Section ti, is essential for 
estimating the precision with which cosmological parameters 
can be determined from the CMB. 

Intermediate angle experiments have detected temper- 
ature anisotropies which are consistent with a primordial 
origin and in rough agreement with adiabatic theory pre- 
dictions. However, the accuracy and robustness of the re- 
sults does not yet allow strong conclusions to be drawn, 
even when experimental results are combined together (e.g.. 
Bond 1996, Bond & Jafi'e 1997, Lineweaver et al. 1997, 
Rocha & Hancock 1997). An experiment with an angular 
resolution of Of^hm ~ 5' can yield useful information about 
the CMB spectrumrj Ci up to multipoles beyond a Gaussian 
flhering scale Is = V8 In 26 J^^^ ~ 2000. If the sky cover- 
age is complete, each multipole is statistically independent. 
However, it is clear from visual inspection {e.g. Figure hn) 



* The temperature power spectrum is defined as the expectation 
value Cf = (|a£mP)i where the coefficients a^„ are defined by a 
spherical harmonic expansion of the temperature anisotropies on 
the celestial sphere AT/T = Y^ o,lmYim{9, <t>)- 
T h is the present value of the Hubble parameter Hq in units of 
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that a typical inflationary Ce curve is smooth and can be 
specified accurately by many fewer than 2000 parameters. 
It is therefore not obvious a priori to what extent the typi- 
cally 10-15 parameters specifying an inflationary model can 
be disentangled by a particular set of measurements. Evi- 
dently, a detailed calculation is necessary as has been per- 
formed in an important paper by Jungman et al. (1996). 
However, as described in the next section, an accurate as- 
sessment of the degeneracies between cosmological parame- 
ters imposed by high resolution CMB experiments requires a 
precise numerical calculation, rather than the semi-analytic 
approach adopted by Jungman et al. 



2 PARAMETER ESTIMATION WITH PRIORS 

2.1 The covariance matrix 

Errors on a set of cosmological parameters s = {si},i = 
1, . . . ,n are estimated using Bayes theorem, which updates 
the prior probability P(s|prior) for the parameters with a 
likelihood function jC{s): P{s) oc £(s)P(s|prior). (Although 
a uniform prior seems to be the least prejudicial, there are 
certain non-debatable restrictions and plausible constraints 
which are reasonabl e to include as prior information, as 
discussed in Section 
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If the errors <5s = s — So about 
the mean So = (s) are small, then an expansion of ln£ to 
quadratic order about the maximum gives:rl 
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adopting a Gaussian approximation for the beam function 
Bel- In this approximation, with a uniform prior the covari- 
ance matrix M = {5s5s') is the inverse of the Fisher infor- 
mation matrix F (e.g., Tegmark et al. 1997), the 1-sigma 

1/2 

error on Si is ai = M^^ , and the correlation coefficient be- 
tween parameters i and j is rij = Mij/{ai(jj). Equation (El) 
gives the standard error AC^ on the estimate of C^ for an 
experiment with A'^ frequency channels c, angular resolu- 
tion 9c,fmhm, and sensitivity ac,pix per resolution element 
{Ocj-whm X dcjwhm pixcl) , which Samples a fraction /aky of 
the sky. 



100 km s~^Mpc~^. The parameters Qi denote the cosmological 
densities of various components (defined in Section 2.2) in units 
of the critical density. 

•f The likelihood also has terms which depend upon the devia- 
tions of the specific realization of the theory represented by the 
observations from the ensemble- averaged values used here (Bond 
1996). We have tested these effects in detail on parameter esti- 
mation and find them to be small and fully consistent with the 
error estimates we quote in the Tables. 
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Figure 1. Temperature power spectra for a standard CDM 
(SCDM) model with h = 0.5, Om = 1 and Ca = and an open 
CDM (OCDM) model with h = 0.6 , fim = 0.33 and Qa = 0. 
The points show a quadratic power spectrum estimation of Ce 
(in bands of width 5% in £), along with the one sigma error, from 
simulated CMB skies observed by the MAP and Planck satellites 
described in Section 3 and Table 1. The errors in the left panels 
use the MAP"*" parameters for the three highest frequency chan- 
nels and two years of observing, while those in the right panels 
use the four lowest frequency Planck HFI channels. Planck can 
follow the theoretical curves quite precisely far down the damp- 
ing tail for both models. In the bottom right panel, a C^ curve 
with Qrn = 0.15, Ca = 0.44 and h = 0.9 almost degenerate with 
the OCDM model is superposed. The inset, showing these two 
and also a Qm = 0.10, Qa = 0.58 and h = 1.1 model, demon- 
strate that the angle-distance degeneracy relation (see text) used 
to define these models is almost exact; these models can be dis- 
tinguished only at low multipoles no matter how precise the CMB 
experiment. 



A key simplifying assumption in deriving equation (0) 
is that the noise is homogeneous and isotropic. Provided 
that the sampling varies slowly, accurate error estimates can 
be obtained using the average weight for the experiment. 
Equation (bl) is exact if fg^y = 1 (e.g., Knox 1995) and is 
approximately correct at multipoles £ > Icut corresponding 
to angular scales small compared to the dimensions 2TT/£cut 
of an incomplete sky map. 

The accuracy of cosmological parameter estimation via 
the covariance matrix approach depends on: (1) the validity 
of the Gaussian approximation to the likelihood function; 
(2) the number and choice of the parameters s defining the 
theoretical model; (3) the parameters So of the target model; 
(4) the numerical accuracy of the derivatives of Ce; (5) the 
inclusion of prior constraints on the parameters s; (6) sys- 
tematic errors in estimates of Ce caused by Galactic and 
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extragalactic foregrounds. For high resolution experiments 
which tightly constrain many of the cosmological parame- 
ters, a Gaussian approximation about the maximum like- 
lihood should be quite good (Knox 1995, Spergel private 
communication), although positivity and other constraints 
can truncate general excursions in the likelihood space. We 
ignore systematic errors caused by foreground subtraction, 
since over much of the sky these are very likely to be much 
smaller than the variance of equation (bh (see e.g. Tegmark 
and Efstathiou 1996 for a discussion of foreground removal 
from CMB maps). Here we consider the remaining four 
points. 



2.2 Choice of variables 

Parameters describing the theoretical angular power spectra 
include those for initial conditions and those characterizing 
the transport of radiation through photon decoupling to the 
present. If we were to allow all possible variations, the count 
of parameters could easily exceed 20; in our analysis, we 
use < 11 variables. We characterize the initial fluctuation 
spectra by an amplitude and a spectral index (tilt) for the 



1/2, 



^1/2, 



itational potential fluctuations, V.J {k„ 
wave fluctuations, VQ^^k. 



scalar and tensor components, "P^' (fc„) and Ws, Vj^^i^") 
and rit. The primordial amplitude parameters for the grav- 

and the gravity 

Gwi^")' ^^^ chosen here to be normal- 
ized at a wavenumber corresponding to the horizon scale. 
The fluctuations arising from inflation could be much more 
complicated, requiring, for example, the parameterization of 
variations of the spectral indices with wavenumber k, inclu- 
sion of isocurvature as well as adiabatic components in the 
scalar perturbations, and possibly of non-Gaussian features. 

At the time of decoupling, the key parameters deter- 
mining the temperature power spectrum are the densities of 
various types of matter, the expansion rate, the sound speed, 
and the damping rate; all of these depend only on the density 
parameters uij = Qjh^, where j — b,cdm,hdm,'y,eru refers 
to baryons, cold dark matter, hot dark matter, and the var- 
ious relativistic particles present then, such as photons and 
relativistic neutrinos. The Hubble parameter at that time 
only depends upon cOm — uJb + ^cdm + i-^hdm. (if the massive 
neutrinos were nonrelativistic thenPI) and ujsr — lo^ + uj^rv. 

The transport to an angular structure of scale (T^ now 
from the post-decoupling spatial pattern of temperature 
fluctuations of comoving scale k~^ depends on the cosmo- 
logical angle-distance relation, t ~ fcT?., where 
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for an open universe (e.g.. Bond & Efstathiou 1984). Here 
OJA = f^Ah^ parameterizes the energy density associated 
with a cosmological constant A (Qa = A/(3//o)) and LOh = 
(1 — f2o)h^ = (1 — flrn — fiA)h^ parameterizes the energy 
associated with the mean curvature of the universe. This re- 
sults in a degeneracy along (5(ajm 71) — lines, which leads 



to a linear relation between StOk and Scja for fixed cOm, with 
coefficients that depend upon the explicit target model. The 
angular pattern we observe also depends upon the change 
of the gravitational metric in time between post-decoupling 
and the present, which breaks this degeneracy. However, this 
late-time integrated Sachs- Wolfe effect influences only low 
multipoles which have a large cosmic variance. Thus, there 
exists one combination of variables which cannot be deter- 
mined accurately from CMB observations alone, even with 
a high precision experiment such as the Planck Surveyor, as 
the lower right panel of Fig. hi illustrates. 

Some parameters are tightly constrained by measure- 



ments other than CMB anisotropies. For example. 



de- 



pends on the temperature To of the CMB, ujeri' depends as 
well on the number of massless neutrino types; the Cf 's also 
depend on the helium abundance, parameterized by Yhe. 
Rather than allow such parameters complete freedom, we 
use the prior probabilities to restrict their allowed varia- 
tions. (Since the experimental errors on Yne, To and A^^ are 
small PI, they have a weak effect on other cosmological pa- 
rameters and hence we include only Yhe in our analysis to 
illustrate the methodology.) 

There also could be many parameters needed to charac- 
terize the ionization history of the Universe; here we use the 
Compton optical depth tc from the present to the redshift 
of reheating, assuming full ionization. We therefore analyse 
a maximum of 11 parameters in this paper: Yne, tc, 4 initial 
condition parameters and 5 density parameters uuj. 

For a given model, the amplitudes of the scalar and ten- 
sor power spectra are uniquely related to the observed am- 
plitude of the CMB power spectrum and that of the present 
day mass fluctuations (characterised, for example, by the 
rms density fluctuation in spheres of radius 8h~^ Mpc, as). 

1/2 

For example, Jungman et al. used the quadrupole Cj 
to flx the amplitude of the fluctuation spectrum. We use 
{£{£ + l)Ci/{2-jT))g , an average over the total band B of 
multipoles that is accessible to the experiment, since this 
is most accurately determined. However the normalization 
parameters ag and pj {k„) are of sufficient interest that 
we also show the accuracies with which these can be de- 
termined. To characterize the tensor amplitude, we use 



rt 



C 



(T) ,r<(.S) 



^1/2, 



1/2/ 



instead of Vj^^ikn) 



In inflation models, 

^Gw('^")/^<i:/^(^") ^'^ simply related to the tensor tilt, with 
small corrections dependent upon the scalar and tensor tilts, 
so one of the four initial fluctuation parameters is a func- 
tion of the other three, here chosen to be nt- rts also depends 
upon other cosmological parameters as well as the tilts (e.g.. 
Bond 1996, equation (6.38)). 

Any parameter set which defines a coordinate system on 
the likelihood surface is a viable set. However, parameters for 
which the Fisher matrix analysis is particularly well suited 
are those for which the first order expansion Ce = Cf (so) -f 
{dCt{so)/ds) • (s — So) is more accurate than the sampling 
variance AC^ for parameters s that he within a few standard 
deviations from the target set So . The set of variables that we 
have adopted gives acceptably high accuracy for the CMB 
experiments described in Section 3. 



S Massive neutrinos become nonrelativistic below a redshift ~ 
1700(mi//ev), where iOhdm ~ 0.01(m^/ev)Afmi^, Nmv is the num- 
ber of neutrino species of mass m^,; for small my, the neutrinos 
may be relativistic or semirelativistic at decoupling. 



1 Yue. = 0.23 ± 0.01 (Pagel et al. 1992), To = 2.728 ± 0.004 
(Fixsen et al. 1997), N^ = 2.991±0.016 (LEP Electroweak Work- 
ing Group 1995). 
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2.3 Choice of target models 

We analyze three spatially flat {flk = 0) target models and 
one with negative curvature. For the canonical flat universe 
we use a standard CDM model (SCDM) with the following 
parameters: ris = 1, n* = 0, Qm = 1.0, fit = 0.05, rts = 0, 
ficdm = 0.95, nhdm = 0, h = 0.5, Tc = 0, Yh, = 0.23; 
it has erg — 1.2, normalized to match COBE DMR, but 
with too many clusters to match the observations. The open 
model has h = 0.6, flm = 0.33, Qb ~ 0.035, Qcdm = 0.30 
and (Tg = 0.44, the COBE-normalization, which has too 
few clusters. We also discuss results for two other DMR- 
normalized spatially flat models that more closely match 
observations: an HCDM model with 2 species of massive 
neutrinos, ilhdm ~ 0.2, m,^ = 2.4 eV, h = 0.5 and erg = 0.77; 
a ACDM model with ^a = 0.67, h = 0.7 and ag = 1.1. (AU 
models have a 13 Gyr cosmological age.) 



2.4 Accuracy of the power spectrum derivatives 

Computational errors in the derivatives of Ce can lead to 
large errors in the covariance matrix. We distinguish be- 
tween two classes of error, one caused by inadequate semi- 
analytic approximations to the Cf and the second caused 
by numerical errors in Ce and its derivatives computed 
from linear Boltzmann transport codes. The Ce accuracy 
must be 1% or better, especially for high resolution exper- 
iments probing multipoles £ > 1000 where the expected 
random errors on each individual multipole become < 3%. 
Errors which are weakly correlated with physical parame- 
ters are particularly serious since these can artificially break 
real near-degeneracies between cosmological parameters and 
lead to overly optimistic error estimates, sometimes by an or- 
der of magnitude or more. Extreme care is therefore required 
in computing the Ce derivatives. For this work we have used 
derivatives calculated with two Boltzmann transport codes, 
an updated version of the multipole code described by Bond 
and Efstathiou (1987) generalized to low density universes 
and including tensor components (Bond 1996 and references 
therein), and the fast path- history code developed recently 
by Seljak and Zaldarriaga (1996). Generally the C^'s from 
these codes agree to better than 1%. We use intervals of typ- 
ically 1 — 5% in the parameters s in computing numerical 
derivatives of C^, i.e. small enough that the derivatives are 
insensitive to the size of the interval, but large enough that 
they are unaffected by numerical errors in the Ce coefficients. 
The primary limitation on the error estimates for Table 
should be the Ce linearization assumption made in deriving 
equation |l|, although we believe that better than percent 
level accuracy in Ce is needed to achieve high precision in 
the nearly degenerate directions of parameter space. The dif- 
ferences between our error estimates and analogous results 
of Jungman et al. , which are large for some parameters, are 
caused primarily by their use of semi-analytic approxima- 
tions to calculate Ct and its derivatives. Zaldarriaga et al. 
1997 have undertaken a similar analysis and come to sim- 
ilar conclusions as those presented here. They also showed 
that polarization information can improve the accuracy of 
some variables, e.g., tc, if foregrounds are ignored. Little is 
known about how the polarization of foregrounds will com- 
promise the relatively weak polarization signal of primary 
anisotropies, especially at low £ where much of the improve- 



ment comes from. See also related work on forecasting errors 
by Mageuijo and Hobson (1997). 



2.5 Inclusion of prior information on parameters 

We have mentioned that a non-uniform P(s|prior) is par- 
ticularly useful for parameters such as Yhb and To, where 
other experiments restrict their values to much higher ac- 
curacy than can be achieved from CMB experiments alone. 
For other parameters, e.g. ui^ and uOcdm, it may be that the 
distribution derived from a CMB experiment is much nar- 
rower than any reasonable prior distribution, in which case 
we gain little by including prior information. There are also 
intermediate cases where prior information can help break 
degeneracies between parameters estimated from the CMB 
alone. We approximate the prior distribution of parameter 
values by a Gaussian distribution with covariance matrix T, 
so the covariance matrix of parameter values including prior 
information is given by 



M = (F + T"^)"^ ff P(s|prior) oc e 



-itfstT-lis 



(4) 



{e.g. Knox 1995), where F is the Fisher matrix (hi). If we are 
interested in the error bars on Si irrespective of the values of 
the other variables, we would marginalize over these, with 

1/2 

error ai — A/./ for the Gaussian case. For most of the 
entries in Table H we use no prior at all ('noP'), except 
for Yhb where indicated. When priors are used, we adopted 
a diagonal covariance matrix Tij with the following values 



\l/2, 



,1/2 



for vTii: 0.3 on the normalization 5{Ce)g /{Ce)g^, 0.5 on 
Us, 2 on rts, 0.075 on cub, 1 on lj,„, 1 on uja, 1 on ujk, 0.5 on 
^hdm , and 1 on re . Some variables are restricted for physical 
reasons to lie within a certain range e.g., tc and rta must be 
positive. Such constraints can be incorporated into the prior, 
but at the expense of more complicated expressions after 
marginalization over these constrained variables. In some 
cases, imposing physical restrictions can lead to a factor of 
two or more improvement in the accuracy of the parameter 
estimates. 

Generally the errors in the parameters will be correlated 
through nondiagonal components of (F -I- T~^)~^. Linear 
combinations of the parameters which are uncorrelated can 
be found by diagonalizing (F-|-T~^). When the eigenvalues 
of (F -I- T~^) are rank ordered, from highest to lowest, the 
variable combinations corresponding to high values will be 
very accurately determined, while those for the lowest may 
be very poorly determined, representing the most degenerate 
directions in parameter space. In Tables 2 and 3 below we list 
the number of parameter combinations that are determined 
within a ±0.01 and ±0.1 accuracy. 



3 COSMOLOGICAL PARAMETER ERRORS 
FROM MAP AND PLANCK 

3.1 The CMB poAver spectrum estimated from 
MAP and Planck 

In this section we apply the above machinery to deter- 
mine the accuracy of cosmological parameter estimation 
from two satellite experiments: the MAP satellite selected by 
NASA (Bennett et al. 1996b) and the Planck Surveyor Mis- 
sion (formerly named COBRAS/SAMBA) selected by ESA 
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(Bersanelli et al. 1996). These satellites offer examples of the 
best that is likely to be achieved in the next decade. Ground 
based and balloon borne experiments will certainly continue 
to provide improved constraints on cosmological parameters 
over this timescale, and so we also analyze a sample long 
duration balloon experiment (LDB). 

The specifications adopted for MAP and Planck are 
listed in Table 1 and have been computed from the infor- 
mation provided on the respective WWW pages for the two 
missions. Although indicative of the expected performance 
of each satellite at the time of writing, these are likely to 
evolve. Of the 5 HEMT channels for MAP, we assume that 
the 3 highest frequency channels, at 40, 60 and 90 GHz, will 
be dominated by the primary cosmological signal. We also 
present the gains that result from a 25% improvement in 
angular resolution at all frequencies and 2 years of observ- 
ing time (we denote these specifications by MAP+). Such 
an improvement is now expected for MAP (Page, private 
communication). Planck will have two detector arrays, a 
Low Frequency Instrument (LFI) using HEMTs and a High 
Frequency Instrument (HFI) using bolometers. The current 
design of the HFI incorporates an additional channel at 
100 GHz in addition to channels at 150, 217 and 353 GHz; we 
have adopted parameters as listed in Table 2 for these four 
channels. We also present results for the 3 highest resolution 
channels in the current design of the Planck LFI which has 
an expected performance that is significantly improved over 
those given by Bersanelli et al. (1996). 

For each multipole £, the computational procedure auto- 
matically rotates the channels into a linear combination op- 
timal for the CMB. In practice, a more sophisticated treat- 
ment would be required in practice to remove Galactic and 
extragalactic foregrounds. It is beyond the scope of this let- 
ter to assess the systematic errors in parameter estimates 
arising from inaccurate foreground subtraction. We there- 
fore simply assume that Galactic foregrounds are negligible 
over a fraction of the sky fsky = 0.65, similar to the 'clean' 
sky area adopted in most analyses of the COBE power spec- 
trum. 

Figure hi shows examples of Ce estimates from one re- 
alization of the SCDM and OCDM target models. In this 
figure, the estimated power spectra have been averaged over 
5% wide bands in £. At the resolution of MAP there is very 
little useful information beyond i ~ 800 (the third acoustic 
peak for the spatially fiat models), whereas Planck samples 
the power-spectrum at close to the theoretical variance limit 
to multipoles £ ~ 2500. The consequences of these differ- 
ences for CMB parameter estimation are described in the 
next section. 



3.2 Accuracy of cosmological parameters 

Results of the analysis for the sample LDB experiment and 
for the MAP and Planck satellites are given in Table Q. For 
the LDB example we adopt specifications for TopHat (Meyer 
et al. 1997), which would cover 4.3% of the sky with an error 
of ISfjK per 20' pixel (which includes an allowance for the 
extra error incurred in removing foregrounds). We assume 
65% of this area will be usable. Because the sky coverage is 
so limited, COBE's DMR is added to improve the baseline 
in i covered and thus the accuracy of parameter estimation. 
(When we apply our analysis to DMR alone, using the av- 



erage noise in the 53-1-90-1-31 GHz map, apix ~ 30fj,K per 
5.2° pixel, we predict the bandpower would be determined 
to 0.07 and n, to 0.20 for SCDM, 0.07 and 0.28 for OCDM, 
if only these two parameters are used, in agreement with 
what is actually found (Section 1); with all 9 parameters, the 
errors grow, but the first and second eigenparameter com- 
binations have error estimates of 0.09 and 0.20 for SCDM.) 
When treating two experiments at very different scales, such 
as DMR and the long duration balloon experiment example 
shown, the log likelihoods (and Fisher matrices) just add. 
Instead of TopHat, we could have chosen any of the other 
bolometer-based LDBs, such as Boomerang (Lange et al. 
1997) or MAXIMA (Richards et al. 1997) or even HEMT- 
based LDBs, such as BEAST (Lubin et al. 1997) and derived 
similar error forecasts. 

For most parameters the inclusion of prior constraints 
on their variation have no effect, particularly for a high pre- 
cision experiment like Planck. Even for MAP the inclusion 
of priors has little impact, except for variables such as Yh^ 
which are poorly constrained from the experiment alone. If 
the helium abundance is allowed to float freely, it has a sub- 
stantial effect on the other parameters; however, limiting its 
value to be 0.23 ± 0.02 results in little impact on the other 
numbers. For the LDB-I-DMR case, the errors are susbstan- 
tially larger with no controlling priors. 

As expected, the parameters have correlations among 
themselves that range from weak to very strong in all mod- 
els, and can differ from experiment to experiment as well as 
model to model. The power amplitude {Ce)g and tc have 
a correlation coefficient about 90% for SCDM for Planck 
and MAP. The most highly correlated are cuk and loa , as ex- 
pected from the angle-distance near-degeneracy. In the Ta- 
bles, the QAii-^ numbers are determined with Q^h^ fixed, 
and the fifch^ numbers are determined with f^Ah^ fixed; 
the other parameters are relatively insensitive to fixing ei- 
ther, or neither. Thus, although our estimates of errors after 
marginalization are gratifyingly small for many parameters, 
especially for the specifications of Planck, they are large in 
other cases {e.g. 5f2Ah^). Error estimates in square brack- 
ets are those obtained when the most correlated component 
for that variable is constrained to be the target value. A 
more natural way to deal with strong correlations between 
variables is to perform a principle component analysis in 
parameter space, rank-ordering linear combinations of pa- 
rameters, as described in Section |2.5| : some linear combina- 
tions are determined exquisitely well and some are less well 
determined because of near-degeneracies as is illustrated in 
the Tables. The Tables also show values obtained in round 
brackets when positivity constraints on parameters such as 
re are used. 

The 5h/h shown are determined from h^ — y^ . loj , 
hence it is a derived rather than fundamental quantity. How- 
ever, h errors depend upon what is kept fixed and what is 
varied. Thus we can use h to replace one of OJ™, oja, cJk, with 
the other two to be marginalized. In that case, the error es- 
timate would be (5h/h = 0.5Sujj/h^. 

We find that the estimated errors on parameters are 
sensitive to their input target-model values. Table H shows 
results for two other flk = models, a ACDM model and 
an HCDM model. This illustrates the sensitivity of param- 
eter error estimation to relatively modest changes in the 
target Ce. In interpreting these tables it is also important 
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to take into account the restrictions that we have imposed 
on the models. The OCDM model error estimates are de- 
rived a ssuming there is no tensor component , Including it 



Bersanelli, M. et al. 1996, 

The Phase A Study for an ES A M3 
D/SCII'96)3: Planck home page, ptto: //astro. cstcc.csa.nl/SA- 



COBRAS/SAMBA, 
Mission. ESA Report 



has Irt.tlp pffppt nn t.Vip rptjiilt.t;- pvpn thp mnst. rnrrplatpH thp 



gcncral/Projccts/Cobras/cobras.htnij 



amplitude, tc and Us, are only slightly affected. The angle- 
distance scaling ensures that the tensor power spectrum does 
not fall off until higher i than in the Qk = cases, and this 
leads to a substantial improvement in Srts- 

We have also found that the errors on nt and rta are 
extremely dependent on the input rts if they are allowed 
to vary independently (Knox and Turner 1994. Knox 1995. 



Efstath ,ou 1997). However, h and the various matter den- 
sities, ujcdm etc., are insensitive to the tensor spectrum for 
reasonable values of rta i$ 2. In open universes, features in 
the power-spectrum shift to larger multipoles according to 
the angle- distance relation, roughly as C{£) ~* C{£/Q,q ); 
thus, for low fio, high resolution is required to determine pa- 
rameters which affect the Doppler peak structure {e.g. h and 
the various w). The relative accuracies of the parameters are 
less sensitive to variations of Qb and h. 



4 CONCLUSIONS 

In summary, we have described how to compute the errors 
in the estimation of cosmological parameters from measure- 
ments of the CMB power spectrum at a number of frequen- 
cies with different angular resolutions and sensitivities. We 
have also shown how prior information on the values of pa- 
rameters can be incorporated into the analysis and described 
some of the pitfalls of this type of analysis that can arise if 
inaccurate derivatives of the C^'s are used and if poor pa- 
rameter choices are adopted. 

We have applied our machinery to the MAP and Planck 
satellites and find that these missions are capable of deter- 
mining fundamental cosmological parameters to an accuracy 
that far exceeds that from conventional astronomical tech- 
niques. In particular, Planck is capable of determining the 
Hubble constant and the baryon density parameter fit, to a 
precision of a few percent or better for each of the target 
models listed in Tables H and H. However, some parameter 
combinations are poorly determined by CMB observations 
alone as described in Section E^ and Section & Nevertheless, 
despite this caveat, it is evident from this work that accu- 
rate CMB observations have the potential to revolutionize 
our knowledge of the key cosmological parameters describing 
our Universe. 
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Table 1. Experimental parameters adopted for this study, u,./^ denotes the central frequency 
of the channel, Of^jhrn the resolution, cTpin,. the pixel sensitivity (in AT/T) per S?„^^ res- 
olution element, £s is the Gaussian beam filter scale, and ui~^ is the noise power for each 
channel. 



MAP (first 


3 used) 












!/,h(GHz) 


90 


60 


40 


(30) 


(22) 




9fwhm 


18' 


23' 


32' 


(39') 


(54') 




fypix/w-^ 


13 


9.9 


7.3 


(6) 


(4) 




^-1/10-15 


4.5 


4.5 


4.5 








is 


465 


345 


255 








MAP+ (w- 


1 xO.5 


(2 yrs), 


Ofwhm X 


0.75, o-p 


^x X 0.94) 




w-i/lO"^^ 


2.3 


2.3 


2.3 








Is 


620 


460 


340 








Planck HFI 


(first 4 


used) 










^,h(GHz) 


100 


150 


220 


350 


(545) 


(857) 


Ofwhm 


14.5' 


10' 


6.6' 


4.1' 


(4.4') 


(4.4') 


<Tpi./10-6 


1.3 


1.3 


1.2 


16 


(77) 


(4166) 


w-i/io-is 


0.028 


0.015 


0.005 


0.35 






is 


560 


800 


1225 


1970 






Planck LFI 


(first 3 


used) 










i^ch (GHz) 


100 


65 


44 


(30) 






9fwhm 


10' 


16' 


23' 


(34') 






^pix /lO-^' 


6.2 


3.7 


2.6 


(1.8) 






w-ViO"^'' 


0.33 


0.29 


0.29 








Is 


810 


505 


350 









© 0000 RAS, MNRAS 000, 



Bond, Efstathiou & Tegmark 



Table 2. Parameter estimation for the 2 models shown. The standard CDM model has 
fim = 1, h = 0.5. The open CDM model has Qm = 0.33, h = 0.6. The columns refer to 
experimental parameters for a fiducial Long Duration Balloon Experiment (LDB), and MAP 
and Planck, with parameters listed in Table 1. MAP"^ assumes a 25% improvement in beams 
and 2 years as opposed to 1 year of observing. Hyvh^ is determined with f2j;h^ fixed with 
a prior, and Hj-h^ is determined with Hyvh^ fixed; most other parameters (except for uj^n) 
are insensitive to fixing either, or neither. Values in square brackets indicate the reduced 
error when the dominant correlated variable is fixed. Values in circular brackets indicate 
what happens when a positivity constraint is imposed. Only selected cases with significant 
variations are shown. In the satellite cases, £cut = 2 was chosen; £cut = 12 was used for the 
LDB experiment, with £cut = 3 for DMR which was analyzed with it. The LDB parameters 
are based upon observing for ten days with the TopHat experiment, and assuming 65% of 
a 24° radius patch will be usable. Because of the limited sky coverage the LDB likelihood 
is combined with the DMR likelihood to constrain low ^'s. DMR has jU~^/10~ = 950. 
{w is the total weight, eq.(3).) Priors are important for the LDB+DMR column, make small 
differences in the MAP column, and essentially none in the rest. P means the prior constraint 
controls that parameter's value. Sample values for Helium with no prior are also shown in 
curly brackets. 



Par am 


LDB 


MAP 


MAP+ 


Planck(LFI) 


Planck(HFI) 




P 


noP 


noP 


noP 


noP 


«j-Vio~^^ 


1.5 


1.5 


.77 


0.10 


.0033 


fsky 


.028 


.65 


.65 


.65 


.65 


SCDM MODEL 


Orthogonal Parameter Combinations within e 


e < 0.01 


1/9 


2/9 


3/9 


3/9 


5/9 


e<0.1 


4/9 


6/9 


6/9 


6/9 


7/9 


Single Parameter Errors from Marginalizing Others 


5{C,fJ^/{CifJ^ 


.022 


.019 (.012) 


.017 


.019 


.015 [.007] 


SUs 


.13 


.06 (.03) 


.04 


.01 


.006 


Srts 


.89 


.38 (.30) 


.24 


.13 


.09 


snth^/n^ho^ 


.23 


.09 (.06) 


.05 


.016 


.006 


sn^h^/ho^ 


.33 


.18 (.11) 


.10 


.04 


.02 


<5f7Ah2/ho" 


.84 


.49 (.35) 


.28 


.14 


.05 


•^f^hdmh^ /ho ^ 


.25P 


.07 


.05 


.04 


.02 


re 


.30 


.22 


.19 


.18 


.16 


SYHe/YHe 


.09P 


.09P {1.4} 


.09P {.59} 


.08P {.19} 


.07P {.10} 


Sag /erg 


.28 


.28 


.23 


.21 


.18 [.06] 


&vy\k^)/vi'^ 


(kn) .24 


.24 


.19 


.17 


.15 [.02] 


5f7feh2/h0^ 


.14 


.07 


.04 


.02 


.007 


<5h/h 


.33 


.19 


.11 


.06 


.02 


OPEN CDM MODEL 


Orthogonal Parameter Combinations within e 


e<0.01 
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2/7 


2/7 


3/7 
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e < 0.1 
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Single Parameter Errors from Marginalizing Others 
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.02 [.016] 


.02 


.02 
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.01 


.003 
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.70 
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.03 


.008 
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.41 
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.006 


^QAhVho" 


1.2 


.31 
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T"C 


.24 
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.07 
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<5Qfeh2/ho" 
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.10 


.07 


.03 


.005 
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.76 


.26 


.18 


.07 


.013 
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Table 3. Parameter estimation for the untilted ACDM model, with Q/^ = 0.66, Q^ = 0, h = 0.70, 
and the untilted HCDM model, with h = 0.50, 2 species of neutrinos with the same mass giving 
^hdm = 0-2, Ha = and Qj. = 0. Errors on Cl/^h^ and O^h^ are independently determined, as 
described in the text. The columns are as in Table 2. 



Par am 


LDB MAP MAP+ 
P noP noP 


COSA(LFI) 
noP 


COSA(HFI) 
noP 


Jsky 


1.5 1.5 .77 
.028 .65 .65 


0.10 
.65 


.0033 
.65 


ACDM MODEL 


Orthogonal Parameter Combinations within e 


£ < 0.01 
£ <0.1 
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Single Parameter Errors from Marginalizing Others 
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